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Abstract 

The m-T-dependencies of the femto-scales, the so-called interferometry and source radii, are 
investigated within the hydrokinetic model for different types of particle pairs — pion-pion, kaon- 
kaon, proton-proton and proton-lambda, — produced in Pb+Pb and p+p collisions at the LHC. In 
particular, such property of the femto-scales momentum behavior as my-scaling is studied for the 
systems with (w) and without (w/o) intensive transverse flow, and also w and w/o re-scattering at 
the final afterburner stage of the matter evolution. The detail spatiotemporal description obtained 
within hydrokinetic model is compared with the simple analytical results for the spectra and 
longitudinal interferometry radii depending on the effective temperature on the hypersurface of 
maximal emission, proper time of such emission, and intensity of transverse flow. The derivation 
of the corresponding analytical formulas and discussion about a possibility for their utilization by 
the experimentalists for the simple femtoscopy data analysis is the main aim of this theoretical 
investigation. 
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I. INTRODUCTION 


The spatiotemporal structures of particle emission in nucleus-nucleus, proton-(anti)pro- 
ton and proton-nucleus collisions are essentially defined by the dynamics of the collision 
processes 1-4]. Therefore the correlation interferometry [o|. that measures the femtoscopic 
scales by means of two- (or many-) particle correlations, allows one to study such pro¬ 
cesses experimentally. The corresponding femtoscopic patterns can be presented in dif¬ 
ferent forms. One of them is the ^-momentum dependence of the interferometry radii 
Ri{hr = |Pn + PT 2 I /2), that results from a 3D Gaussian fit in q l = p\ — p 2 to the two- 
particle correlation function C7(q, kr), defined as a ratio of the two-particle spectrum to the 
product of the single-particle ones. The other one is the source function S(r*) 6] reflecting 
the dependence of the pair production on the distance r* between the two emitted particles 
in the rest frame of the pair. Both patterns supplement each other, and a reliable model 
should describe/predict all the mentioned types of the femtoscopic observables, if it contains 
a detailed space-time picture of the collision process. 

It is important to note that the correlation function behavior depends also on the particle 
species. The detailed behavior of this dependence can be used to discriminate between 
different scenarios of matter evolution and particle emission in the collision processes. For 
example, the hydrodynamic picture of A+A collisions for the particular case of negligible 
transverse flow leads to the same m^ 2 behavior 0 of the longitudinal radii Ri{kx ) for the 
pairs of identical pions and kaons, and even leads to the complete m^-scaling in the case of 
common freeze-out S3- „„ 

In Ref. [7] it was found that hydrokinetic model (HKM) 00 predicts strong violation 
of such a scaling between pions and kaons at the LHC, and predicts /cr-scaling (for kx > 
ko ~ 0.4 GeV/c) instead. In this letter we analyze in detail the physical reasons for niT- 
scaling violation. As it turned out that the reasons are quite general, there is a hope that 
the found mj- peculiarities will be confirmed in the LHC experimental comparative analysis 
of pion and kaon femtoscopy data. In addition we predict the m^-behavior for femtoscopic 


1 Here, = m 2 + (|pti + PT 2 I /2) 2 = m 2 + fcf. is the transverse mass of the particle pair. 
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scales in other cases, including the case of meson and baryon source function radii. 

To simplify the theoretical study, we reconsider the analytical results, describing tut (or 
kx) behavior of the femtoscopy scales in a situation with strong transverse flows typical for 
RHIC and LHC energies. We derive a simple analytical formula for such a scenario, which 
fits well the complex HKM results for different hadron pairs. Since the HKM describes 
simultaneously the large number of bulk observables, it gives experimentalists a simple tool 
for estimation of the life-time of expanding fireball, related to the maximal emission of the 
specific hadron pair created in proton or nuclear collisions. 


II. ANALYTICAL MODEL FOR THE INTERFEROMETRY RADII 


The hydrokinetic model (HKM) (see details in Ref. Q) was developed to simulate 
the evolution of matter formed in relativistic heavy-ion collisions and describe/predict bulk 
observables at RHIC and LHC. In Ref. [ll] a good description was reached for pion, kaon, 
(anti)proton and all charged particle spectra at different centralities, as well as for the 
elliptic flows. Concerning the femtoscopic scales, the confirmed HKM prediction 12] about 
a reduction of 4^- ratio in A+A collisions at LHC as compared to RHIC was done; the 

M side - - - 

interferometry radii for pions and kaons at RHIC 
were calculated; the source functions for kaons anc 
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. ll ll ] and LHC {?, H] in A+A collisions 
pions at RHIC were described at the top 
RHIC energy and predicted for the LHC one jl3]. In addition, HKM describes well pion 
interferometry radii in p + p collisions at the LHC energy y/s = 7 TeV if one incorporates 
the quantum uncertainty principle into a quasi-classical event generator [t3 | . 

All the results for A+A collisions are obtained at the same initial conditions at r = 
0.1 fm/c, produced by MC Glauber generator GLISSANDO [150 (with just re-scaled max¬ 
imal initial energy density at the transition from RHIC to LHC), the same type of hydro¬ 
dynamics evolution (however with different baryochemical potentials at RHIC and LHC), 
the same lowest temperature when the local thermal and chemical equilibrium still exist 
in the hydrodynamically expanding medium, T c h = 165 MeV, and with UrQMD cascade 
at the final afterburner stage. Such a self-consistent and reliable simulation of the data, 
including femtoscopy ones, indicates that the spatiotemporal picture of particle emission in 
the model can serve as a reference point for further studies, no matter how much advanced 
future models will be developed. Note, that in HKM and similar models many factors act 


3 















simultaneously, the results can be obtained only by means of time-consuming numerical cal¬ 
culations, and it is not always clear, how each physical factor affects the complete space-time 
picture of the particle emission in the collision process. 

Here we provide some analytical estimates for the spectra and femtoscopy scales, aiming 


to revea’ 
In Refs. 


le main parameters forming the particle emission distribution in space and time. 


. t 

ie 

a 

16 


16] the problem of spectra formation in hydrodynamic approach to A+A collisions 


is considered within the Boltzmann equations. It is shown analytically and illustrated by 
numerical calculations that the particle momentum spectra can be described by the Cooper- 
Frye prescription (CFp) jl3] despite freeze-out is not sharp and has the finite temporal 
width. The latter is equal to the inverse of the particle collision rate at 4-points (t a (r,p), r) 
of the maximal particle emission at a fixed momentum p. The set of these points forms 
the hypersurfaces t a (r,p) which depend on the values of p and typically do not enclose 
completely the initially dense matter. This is an important difference from the standard 
CFp, with a common freeze-out hypersurface (that is, typically, an isotherm) for all the 
momenta p. Also, the well known problem of CFp - negative contributions to the spectra 
from non-space-like parts of the freeze-out hypersurface - is naturally eliminated in this 
improved/generalized prescription: at each momentum p the hypersurface a p is always a 
space-like one. 


In this letter we partially use the results of Refs. 


18 


19] and consider the hypersurface of 


constant Bjorken (proper) time r, a p : r = const, as the hypersurface of maximal emission 


for soft enough particle momenta p. It is shown m that such a hypersurface, limited in 
transverse direction r^, is typical for soft particles. For the particles with hard momenta, 
e. g. Pt > 0.8 GeV/c, the hypersurface of maximal emission is different from r ~ const: 
there is the correlation between the transverse radius and the time of radiation (see details 
,3)- To cut the hypersurface of maximal emission (m. e.) r = r m . e . = const in the 


m 


18 


transverse plane for soft quanta, and at the same time get analytical approximation, we use 


here the Gaussian cutoff factor p(r^) in the way proposed first in Refs, 
the Wigner function for bosons in almost central events has the form: 


18 


19]. Namely, 


fl.eq.(x,p) = 


[exp (f3p • r T ) - ftp) - 1] 1 p(r T ) 


(2tt) 


( 1 ) 


Here and below j3 — 1/T is the inverse of temperature at r m e . = — x 2 L m e , = x-r is 

2D transverse radius-vector, 4-coordinate is x 11 = (r coshp^, r^, Tsinhpi) with longitudinal 
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rapidity tjl = arctanhuL = transverse rapidity rj T = r] T (r T ) = arctanhn T (r T ), 

= (cosh r/i cosh tit, ^ sinh t)t, sinh jjl cosh rjr)- The cutoff 


the hydrodynamic velocity 


{x 

factor p(tt) has the form 

18, 

19 


p( r T ) = exp[—a(cosh? 7 r (r r ) — 1)]. 


( 2 ) 


Here, according to results of 18j, [19] and the interpretation corrected for the improved 


CFp [lfj], the parameter a = R%/R where Rt is the transverse homogeneity length in 
r T (near r T = 0 for small k T ) along the hypersurface r = r m . e . = const, and R v is the 
hydrodynamic length, R v = (^'(ry)) -1 , near the same tt- The small a can be reached at 
very intensive flow when the hydrodynamic length R v is much smaller than the homogeneity 
length Rt (which is, generally speaking, different for different particles), while the upper 
limit is reached in the case of no transverse flow, R v = oo. One can see that the contribu¬ 
tions from particles with high momenta emitted from the hypersurface r = r m . e . = const 
are suppressed by the cutoff factor, since such particles radiate typically from rapidly mov¬ 
ing fluid elements, coshr] T (r T ) 1. Another r( r T ) hypersurface should be used for hard 
particles. 

To get single particle spectra pod 3 N/d 3 p and the correlation function of bosons C(p,q) 
in smoothness & mass shell approximation, we use the generalized Cooper-Frye method as 
it was discussed above: 

d 3 N 


P o- 


d 3 p 


' O'm.e. ( P ) 


da^fi.eqXx^p), 


(3) 


C(p, q) « 1 + 


L m . e .(k) da v kfl fu q .(x, k) exp (iqx) 


(4) 


(fa m .e.(k) da^f Leq .(x, k ) 

where q — Pi — P 2 , k M = [\jm 2 + ( Pl ^ P2 )^, Pl ]( P2 )• Note that in smoothness & mass shell 
approximation k ~ p — {p l -\-p 2 )/2. Then 4-vector q has three independent components, 
which can be selected along the beam axis, qi ong = qi , along the pair transverse momentum 
vector k T , q out = q Q , and along side direction, q side = q s , being orthogonal to both long and 
out directions. 

With all these notations, one can apply the saddle point method (in complex plane) to 
calculate the spectra and correlation function (j4]) in Boltzmann approximation for longitu¬ 


dinally boost-invariant expansion. The results were obtained in Ref. 


20] on the base of the 
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approach using the Wigner function (JT]) |1|], lfj. The remarkable fact is that the behavior 
of the correlation function in long direction depends solely on the parameter a and does not 
depend on the transverse velocity profile at the hypersurface of maximal emission, while, 
as expected, the behavior of the correlation function in out and side directions is model 
dependent and quite sensitive to the details of the velocity profile. Because of this, here we 
analyze and use only the longitudinal projection of the correlation function. Let us introduce 
the value A, associated with the homogeneity length in longitudinal direction in a presence 
of transverse flow A; = t^J ^(1 — v f,) 1 / 2 IS) 12]: 


A 2 = % = —(1 - 4) 1/2 . 
t z m<r 

where the transverse velocity in the saddle point Vt = Jct/(tht + ctT). Here T = T, 


temperature at the hypersurface of maximal emission, r = r m . e .. Then in LCMS 


20 |: 


(5) 


is the 


C(k, q u q s = q 0 = 0) = 1 + 


exp 


2 

1? 


1 - y/l + r 2 X A qj 


Ahr—>• oo 


1 + exp(—A z q t ). (6) 


[1 + T 2 A 4 qf] 3 / 2 

At large m^/T ^ 1 the correlation function has the Gaussian form, the interferometry 
radii coincide with homogeneity lengths in this approximation, that was first obtained for 
pure Bjorken expansion (without transverse flow: vt = 0) in Ref. Q: Ri = A, - r 


In Ref. 21| for the same case (no transverse flow) there was found the correction to the 
interferometry radii at small ratio mx/T « 1, that can be actual for pions: 


Rf — A? -)• X'i 


J<2^) 

KA^y 


(7) 


where K n are modified Bessel functions. As one can see from (|6l) . at not large mr/T ~ 1 the 
correlation function has non-Gaussian form in qi. Correspondingly, at not large m-r/T such a 
radius describes only the peak of the correlation function, where the Gaussian approximation 
in the limit of small qi, C(qi ) = 1 + exp(— Rfqf), can be obtained. The detail analysis of 
different analytical approximations for interferometry radii is done in Ref. 20]. 

The approximation of the correlation function (JS]) at small qi leads to the following analyt¬ 
ical result for the interferometry radii in the case of boost-invariant longitudinal expansion 
and transverse flow with arbitrary velocity profile: 

3 


Rf(k T ) = t 2 A 2 


f + -A 2 


( 8 ) 
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The A is defined by the formula ([5]). The direct comparison with the formula (J7J) for pure 
Bjorken picture (in our case it corresponds to a —> oo) demonstrates only 1 — 3% deviation 
between the two results in all actual kx interval for pions. An advantage of the formula (JH]) is 
that it is derived for the case of transverse flow of any intensity that is especially important 
for LHC energies. 


III. ANALYTIC FITTING OF HKM RESULTS 

The idea of m^-scaling for the interferometry radii of different bosons has been launched 
by the ID hydrodynamic result js]]: Ri cc a/ T/tut■ The formulas (J7J) and (JH]) extended 
for all rriT values in the case of ID expansion (it means a — oo =>- vt — 0 in f-la (THT) ) 
also lead to scaling: Ri is the function of r and ttit/T only. In Fig. |T]we demonstrate the 
corresponding fit to the HKM results for pions and kaons with the standard Cooper-Frye 
prescription at the “freeze-out” hypersurface T = 165 MeV, where we artificially switched 
off the UrQMD re-scattering, leaving only resonance decays. As we see, the best fits are very 
bad: % 2 /ndf = 2562.31 for pions and y 2 /ndf = 585.24 for kaons, and there is no mr-scaling. 
The absence of rriT scaling in full HKM (with re-scattering at the afterburner UrQMD stage) 
for Ri at LHC energies is noticed recently in Ref. [7]. 

One of the reasons for violation of scaling in HKM becomes clear from Fig. [2] where 
the same results of “succise” HKM (without re-scattering) were fitted by the full f-la (J8]) 
that accounts for transverse flow. Here we also use the “freeze-out” temperature 165 MeV. 
The fits for both pions and kaons are quite satisfactory with the same value of temperature 
parameter and times r n = 7.41 frn/c for pions and tk = 7.56 fm/c for kaons. For both 
mesons a = 2.8. These fits illustrate that, as one can see from (JHJ), (J5]), the transverse flow 
destroys scaling, since the interferometry radius depends now on both variables mj/T 
and k T /T. 

The additional reason for m-T-scaling violation could be re-scattering at the afterburner 
stage that also can break the similarity between pion and kaon longitudinal interferometry 
radii. A non-Gaussian shape of the correlation functions also affects the m^-behavior of 
the corresponding femtoscopy scale. To investigate in detail the role of these factors, we 
should analyze the results of full HKM calculation (including the re-scattering stage). In 
order to define the effective T and a parameter values for an analysis of the interferometry 
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FIG. 1. The rriT dependencies of longitudinal femtoscopy radius R[ calculated in HKM model 
without re-scattering stage for K ch K ch and pairs (markers) together with the corresponding 

fits (lines) according to formula © where transverse flow is absent (a = oo). The temperature 
is T = 165 MeV. The results are related to yj snn = 2.76 GeV Pb+Pb collisions at the LHC, 
c = 0 - 5%, \rj\ < 0.8, 0.14 < p T < 1.5 GeV/c. 


radii behavior in full HKM, one can study the corresponding momentum spectra. These 
parameters can be extracted from the combined fit to pion and kaon spectra calculated 
in HKM, accounting for transverse flow. Such a fit can be performed using the formula 
obtained in the same approximation as the result ([6]) [H 


Per 


d 3 N 

d 3 p 


oc exp 


[-( m T /T + a)(l - Vt) 1/2 }. 


(9) 


In Fig. E] one can see the experimental 22] and HKM spectra together with the combined 
pion-kaon fit to HKM points. The fit with a common temperature is performed according 
to (© in p T -range 0.5 — 1.0 GeV/c. The extracted parameter values are the following: 
T = 144 ± 3 MeV, a n = 5.0 ± 3.5 and ax = 2.2 ± 0.7. To maximally reduce the effect of 
the non-Gaussian correlation functions, we restrict the fitting range to q = 0 — 0.04 GeV/c 









FIG. 2. The same as in Fig. Q] but fits account now for the intensive transverse flow, a = 2.8. 
The time parameters for pions and kaons are r n = 7.41 fm/c and tk = 7.56 fm/c correspondingly, 
T n = Tk = 165 MeV. 


when extracting the interferometry radii from full HKM. One can see the m^-dependence 
of the radii extracted in such a way in Fig. [5] together with the analytical fit according to 
Eq. ([HD, where T, a are constrained according to the combined spectra fit. The maximal 
emission times extracted from the fit for pions and kaons are r n = 9.44 ± 0.02 fm/c and 
tk = 12.40 ±0.04 fm/c respectively while the rest of parameters take values T n = 147 MeV, 
T k = 141 MeV, a n = 8.5 and a# = 1.5. For comparison, let us present the pion and kaon 
emission pictures in HKM based on the invariant emission functions G(x,p) = p° d f x fp p ■ In 
Fig. [H we demonstrate the reduced emission function averaged over all momentum angles 


t; pt) = 


p°d 6 N 


dr Q dr s drjdrdpT dy 


r s = 0 


for pions and kaons at 0.2 < px < 0.3 Gev/c in the central 


rapidity bins for both space-time p and momentum y rapidities. Here r 0 is the component of 
along the vector p-r, and r s is its component in the transverse direction orthogonal to p^. 


Such a type of the emission function (in r s and r 0 variables) was used in 18 h 20] to derive 
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FIG. 3. Pion (blue) and kaon (red) momentum spectra in Pb+Pb collisions at the LHC, -y/siviv = 


2.76 GeV. Open squares show experimental values 


221 ]. triangles show HKM results, black lines 


correspond to combined (with the same temperature T) pion and kaon spectra fit according to (|9|). 


the analytical approximation for the correlation function (see Section 2). The fit results for 
the times of maximal emission are in good agreement with the HKM emission picture. One 
can notice that the kaon radiation has two maxima: one happens near r = 10 frn/c and the 
second, broader and less pronounced, is at r ~ 14 — 15 fm/c. This second local maximum 
must be due to decays of A'* (892), having life-time 4 — 5 fm/c, into kaons A' ± (493.7). It leads 
to some kind of mean time for kaons about 12 fm/c. It is worthy noting that the difference 
between this “mean” time and time of pion maximal emission is significantly larger than 
in the case when only resonance decays are taken into account, but not re-scattering. The 
reason is that a free streaming of fast K* with subsequent decay into K ± leads to the 
additional large source, which contribution to the correlation function results mostly in non- 
Gaussian behavior of it. At the same time the re-scattering of K* involves these mesons 
in some kind of collective motion producing thereby the direct link between the time of 
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maximal emission and the value of the longitudinal radii. One can also note that the results 
of fitting the femto-scales extracted from full HKM (accounting for the re-scattering and K* 
decays) as well as combined fitting of the HKM spectra give different T and a parameters 
for pions and kaons, and so one can conclude that in addition to intensive transverse flows 
m T scaling is violated because of the re-scattering and resonance decays. 



r T (fm) 



r T (fm) 


Ho. 007 
- 0.006 
- 0.005 
- 0.004 
0.003 
- 0.002 
^ 0.001 
®0 


FIG. 4. The momentum angle averaged emission functions per units of space-time and momentum 
rapidities g(r,rT,PT) [fm -3 ] (see body text) for pions (a) and kaons (b) obtained from the HKM 
simulations of Pb+Pb collisions at the LHC y/ snn = 2.76 GeV, 0.2 < px < 0.3 GeV/c, \y\ < 0.5, 
c = 0 - 5%. 


Next, if we apply Eq. ([8]) to description of the femtoscopy radii obtained from the fits to 
pion and kaon HKM correlation functions in a broader range (e. g. q = 0 — 0.2 GeV/c) Q, 
covering not only their bell-like parts, but also the flatter ones, we will get unsatisfactory 
model data description for kaons using the fit parameters corresponding to spectra combined 
fit (see Fig. El blue solid line) with T K = 141 MeV, olk = 1-5, and tk = 11.09 ± 0.02 fm/c. 
As for pions, they are described well with parameters restricted by combined spectra fit: 
T n = 141 MeV, a n = 1.82, r n = 10.34 ± 0.06 fm/c. To get more or less adequate fit for the 
kaon longitudinal radius corresponding to wide g-region, one should remove constraints on 
(Xk parameter. The dashed blue line in Fig. El represents such a fit with still constrained 
temperature, T = 144 ± 3 MeV, and free ax and Tx- The resulting fit at T = 146 MeV 
and a = 0.02 provides good data description and gives tk = 12.65 ± 1.58 fm/c, that is in 
good agreement with the previous results for radii extracted at q = 0 — 0.04 GeV/c and the 
emission picture in Fig. [4] Therefore, the combined fit (with common temperature) for the 
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FIG. 5. The same as in Fig. [TJ [2] but radii are calculated in the full HKM model including re¬ 
scattering stage. To reduce the effect of the non-Gaussian correlation functions, we take more 
narrow fitting range for them, q = 0 — 0.04 GeV/c. The fit parameters T and a correspond to 
combined pion and kaon spectra fitting. At T n = 147 MeV, Tx = 141 MeV, a % = 8.5 and ax = 1.5 
extracted maximal emission times are T n = 9.44 ± 0.02 frn/c and tk = 12.40 ± 0.04 fm/c. 

pion and kaon spectra leads to the reliable estimate of the maximal particle emission times 
in the case of fitting only the peak of the non-Gaussian correlation function. The same good 
estimates can be extracted also from the Gaussian fit in large g-diapason, but for kaons ax 
parameter has to be fitted without restrictions coming from common pion and kaon spectra 
fitting. 

The difference as compared to the previous results for “without re-scattering” regime as 
well as an essential difference between the main fit parameters for pions and kaons mean 
that the re-scattering plays an important role in violation of my scaling, in addition to 
transverse flow. The very small a value in kaon fit can indicate that the re-scattering stage 
affects seriously kaon emission picture and the re-scattering contribution to the femtoscopy 
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scale cannot be described basing on the pure hydrodynamic (plus the resonance decays near 
freeze-out hypersurface) approach. The peculiarities of the re-scattering for kaons leads to 
more pronounced non-Gaussian character of the correlation function shape than for pions 
and also to a distortion of the femtoscopic scales dependence. It is found that for 
compensation of the influence of such a factor on the extracted effective time of maximal 
emission, in the case of wide (/-fitting range for kaons, one can eliminate restriction for the 
parameter a coming from combined pion and kaon spectra fit. 



FIG. 6. The same as in Fig. [5] but the radii are extracted from the fits to full HKM correlation 
functions in a wide q range, q = 0 — 0.2 GeV/c. Thus, the significant deviation of the CF shape 
from the Gaussian one leads to a distorted radii m^-behavior. At fitting it can be compensated 
by decreasing ax as compared to that obtained from the spectra fit. Experimental data 
pions are demonstrated for comparison. 


2M for 


In some sense a complementary method of obtaining the information about the space-time 
structure of the system from the correlation measurements, is known as the source imag¬ 


ing 


24 


26] aiming to extract the source function S( r*), that represents the time-integrated 
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distribution of particle emission points separation r* in the pair rest frame. The analysis 
of the experimental data for baryon correlations, often based on the Lednicky-Lyuboshitz 


formula 


27], supposes isotropic Gaussian separation distribution described by single width 


parameter r 0 . This value can be extracted from the Gaussian fit to the angle-averaged 
source function S(r ) = 1/(47r) J^ 77 S(r, 6, 0) sin 9d9d(f>. In this paper we also try to ap¬ 
ply the fitting function (J8]) to pion, kaon, proton and proton-A pairs to describe the nix 
dependence of source radii r 0 extracted from the corresponding angle-averaged source func¬ 
tions S(r*) in HKM by means of Gaussian fits to them, see details in Ref. 13]. The results 
are presented in Fig. [7] At this fitting the temperature T was constrained according to 
the combined pion and kaon spectra fit result T = 144 ± 3 MeV, while a and r were left 
free. The fits are unexpectedly good with the parameters T = 141 MeV for all the par¬ 
ticle pairs, = 11.47 ± 0.03 fm/c, tkk = 11.26 ± 0.04 frn/c, t pp = 11.30 ± 0.13 frn/c, 
T pA — 12.44 ± 0.29 fm/c. The parameter a for all the pairs is about 10 4 and more, i. e. 
actually a —y oo. The latter means that there are no transverse flows in the pair rest frame, 
saddle point in transverse velocity is zero for fluid elements in rest, that is natural. Found 
TOr-scaling (for proton-A pairs rriT = \J( mp + m ± y _[_ with slightly variative multiplier r 

looks very surprising. 


Let us also apply the result ([HI) to p + p collisions at LHC energy y/s = 7 TeV. The 
theoretical results for the pion interferometry radii in p + p collisions with different particle 
multiplicity accounting for uncertainty principle in small systems were presented in Ref. Q. 
The results without re-scattering are almost the same as in full HKM. Therefore we assume 
the temperature of chemical freeze-out to be 165 MeV. In Figs. [8] [9] we fit by f-la (JH]) these 
results obtained in the full HKM model, including re-scattering at the afterburner stage 
and corrections accounting for quantum uncertainty principle. One can see that, as it is 
expected, the values of r in fitting function are much smaller than in Pb+Pb collisions, 
namely, in the multiplicity bin with mean charged particle multiplicity 17.9 we obtained 
r = 1.83 ± 0.02 fm/c, a = 13.9 ± 0.22. At the same time for smaller mean multiplicity 9.2 
both values are smaller: r = 1.60 ± 0.04 fm/c, a = 9.52 ± 0.44. 


14 







9 


71 71 



O L_l_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_L 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

m T (GeV/c) 

FIG. 7. The mx dependencies of tt-k, K ch K ch , pp and pA source radii ro extracted from correspond¬ 
ing angle-averaged source functions calculated in HKM for ^/snn = 2.76 GeV Pb+Pb collisions 
at the LHC, c = 0 — 5%, \p\ < 0.8. Transverse momentum ranges are 0.14 < px < 1.5 GeV/c for 
pions and kaons, 0.7 < px < 4.0 GeV/c for protons and 0.7 < px < 5.0 GeV/c for Lambdas. 

IV. CONCLUSIONS 

The analytical formula for the longitudinal interferometry radii that allows one to extract 
the time of maximal emission for the certain hadron pairs in the case of strong transverse 
flows in A+A and p + p collision processes is proposed. It is compared with detail calcula¬ 
tions and spatiotemporal picture of the particle emission in hydrokinetic model that already 
described well the large number of the bulk observables. Within such an analysis the fac¬ 
tors were clarified affecting the interferometry radii and leading to violation of m^-scaling 
between pions and kaons. They are: strong transverse flows, re-scattering and resonance 
decays at the afterburner stage and non-Gaussian form of the correlation function. It is 
found that for compensation of the influence of the latter factor on the extracted effective 
time of maximal emission for kaons, one may eliminate restriction on the “flow intensity” 
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FIG. 8. The mx dependence of pion longitudinal femtoscopy radius Ri calculated in the full HKM 
model with corrections accounting for quantum uncertainty principle (red markers) together with 
the corresponding fit ([HD (red line). The temperature is T = 165 MeV, thkm = 1.60±0.04 fm/c and 
a = 9.52 ±0.44. The results are related to yfs = 7 TeV pp collisions at the LHC, ( dN c h/dr /) = 9.2. 

parameter a for kaons which follows from combined pion and kaon spectra fit. 

As for the pion interferometry radii in p ± p collisions, the analytical fits to the corre¬ 
sponding HKM calculations give quite reasonable estimates for the times of the maximal 
emission for the collisions with different multiplicities. 

The fitting formula has been applied also to the so-called source radii obtained from the 
Gaussian fit to the source function describing the dependence of the pair production on 
the distance between the two emitted particles in the rest frame of the pair. The pion- 
pion, kaon-kaon, proton-proton and proton-A pairs produced in Pb±Pb collisions at the 
LHC were analyzed within HKM. Very surprisingly, it was found that analytical fit with 
slightly different r for different pairs and with the other common adequate parameters - 
no transverse flow in the pair rest frame and the same temperature taken from combined 
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FIG. 9. The same as in Fig. [8] for collision events with ( dN c h/dr /) = 17.9. The values of the fit 
parameters are thkm = 1-83 ± 0.02 frn/c, a = 13.9 ± 0.22. 

pion and kaon spectra fit — gives very good result and leads to the mj> scaling behavior of 
angle-averaged source radii for all mentioned baryon and meson pairs. 

We hope that the results will be especially interesting for the experimentalists for simple 
fitting of the femtoscopy data and extracting the time of maximal emission in continuously 
emitting fireballs. 
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